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The alternative replica technique which involve summation over all integer momenta of the par- 
tition function and which does not require analytic continuation to non-integer values of the replica 
parameter n is discussed. In terms of this technique (which does not involve any replica symmetry 
' breaking "magic operations") rigorous solution for the average free energy of the Random Energy 

7-H Model is recovered in a very simple way. 
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I. INTRODUCTION 



In recent years there is a renewed interest to the mathematical status of the replica method widely used in disordered 
systems during last four decades (see e.g [HQ)- For the calculation of thermodynamic quantities averaged over disorder 
~v parameters (e.g. average free energy) the method assumes, first, calculation of the averages of an integer n-th power 
of the partition function Z(n), and second, analytic continuation of this function in the replica parameter n from 
I 1 integer to arbitrary non-integer values (and in particular, taking the limit n 0). Usually one is facing difficulties 
at both stages of this program. First of all, in realistic disordered systems the calculations of the replica partition 
^ . function Z(n) can be done only using some kind of approximations, and in this case the status of further analytic 
I t ' continuation in the replica parameter n becomes rather indefinite since the terms neglected at integer n could become 
essential at non-integer n (in particular the limit n — > 0) 0, On the other hand, even in rare cases when the 
derivation of the replica partition function Z(n) can be done exactly, further analytic continuation to non- integer n 
I , appears to be ambiguous. 

The classical example of this situation is provided by the Derrida's Random Energy Model (REM) [5[ . At present 
this is one of the best studied models of spin-glasses (see e.g.@ and references there in) which exhibits non-trivial 
solution. It can be easily shown that in this system the partition function momenta Z(n) growths as exp(n 2 ) at 
large n, and in this case there are exist many different distributions yielding the same values of Z(n), but providing 
different values for the average free energy of the system Q . In this situation the replica solution which is generally 
believed to be correct is obtained via the "magic operations" of the Parisi replica symmetry breaking (RSB) scheme [l( 
(which in the case of REM reduces to the special case of the so called one-step RSB). Unfortunately, this is not more 
\Q • than a heuristic procedure which at present has no rigorous mathematical grounding. On the other hand, it should 
\ be noted that during last decade remarkable progress has been achieved in mathematically rigorous derivations of 
• various results previously obtained in terms of the replica method. A number of rigorous results have been obtained 
| which prove the validity of the cavity method for the entire class of the random satisfiability problems revealing the 
physical phenomena similar to what happens in REM and which are described by the one-step RSB solution (see e.g. 
[2} and references there in). The results obtained in terms of the continuous RSB scheme developed for mean- field 
spin glasses has been also confirmed by independent mathematically rigorous calculations 8], 

Recently a notable progress has been achieved in the replica calculation technique itself 
not require performing analytic continuation to non-integer values of the replica parameter n, and formally makes 
possible to compute an entire free energy distribution function summing over all integer momenta Z(n). In this brief 
communication I would like to present very simple replica calculations which do not involve such tricks like RSB 
"magic operations", and which, nevertheless, recover well known result for the average free energy of REM at all 
temperatures including the phase transition into the low-temperature phase (which is usually called the one-step RSB 
state). 



This technique does 



II. REPLICA TECHNIQUE 



By definition the partition function Z of a given sample is related to its free energy F via 

Z = exp(-/3F) 



(1) 
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The free energy F is denned for a specific realization of the disorder and thus represent a random variable. Taking 
the n-th power of both sides of this relation and performing the disorder averaging we obtain 



Z n = Z(n) = exp(-/3nF) (2) 

where the quantity in the lhs of the above equation is called the replica partition function. The averaging in the rhs of 
the above equation can be represented in terms of the free energy distribution function P(F). In this way we arrive 
to the following general relation between the replica partition function Z(n) and the distribution function of the free 
energy P{F): 

/ + OC 
dFP(F) exp(-/3nF) (3) 
-oo 

The above equation is the bilateral Laplace transform of the function P(F), and it looks as if, at least formally, it 
allows to restore this function via inverse Laplace transform of the replica partition function Z(n): 

r +iao d(nB) 

= i ^El z{n) exp{ p nF) (4) 

In order to do so, first one has to compute Z(n) for an arbitrary integer n and then perform analytical continuation 
of this function from integer to arbitrary complex values of n. This is the standard strategy of the replica method in 
disordered systems where it is well known that very often the procedure of such analytic continuation turns out to be 
rather controversial point [H Q . 

Usually the free energy of a given random system is expected to be the extensive quantity: F = Vf, where V is the 
volume of the system and / is (random) free energy density described by some distribution function Vy(f) (which 
in general depends on the volume V). Substituting this into eq.Q and introducing a new integration parameter 
s = j3nV we get 

V V (J) = [ +lC ° ^Zv(s) exp( s /) (5) 

where Vy(f) = VP(Vf) and Zy(s) = Z(s/j3V). If in the thermodynamic limit V — > oo the phenomenon of 
selfaveraging takes place the limiting free energy distribution function becomes <5-like: limy->oo Pv(f) = S(f — /(/?)) 
where f(j3) is the mean free energy density which is the quantity of the first interest in the disordered systems. 
According to the above equation, this means that the limiting replica partition function is expected to take the form 
limv-j-oo Zy(s) = exp(— s/(/3)) where the parameter s remains finite. Since, by its definition, s = (3nV, this implies 
that in the limit V — >■ oo, the replica parameter must go to zero, n ~ 1/V — > 0. 

The problem is that before taking the thermodynamic limit, the replica partition function Zy{s) has to be computed 
for finite volume V. It is well known that in many cases the finite size distribution functions of random quantities 
are extremely singular objects, and only in the thermodynamic limit they converge to smooth and nice shape. The 
typical example is provided by the eigenvalues distribution functions in the random matrix theory (see e.g.[ll|). For 
that reason it could be easier, instead of the distribution function itself, to study its integral representation, namely, 

/•oo 

W(x) = / dfV(f) (6) 

By definition, the function W(x) gives the probability that the random quantity / is bigger than a given value x. It 
is clear that this function is much more "smooth" object that the distribution function itself: even in the case that 
the finite system size function Vy(f) represents a set delta- functions, its integral representation would be only a kind 
of step-like continuous curve. 

Formally the thermodynamic limit probability function W(x) can be defined as follows: 

W(x) = lim V ^ — f- exp((3nVx) Z n (7) 

n=0 

°° (—i) n 

= lim V ^ — {— exp(/3nVx - /3nVf) 

n=0 



= ^lim exp[— exp(/3F(x — /))] 
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which coincides with the definition, eq.©. Thus, according to eq.([7]), the probability function W(x) can be computed 
in terms of the above replica partition function Z(n) by summing over all replica integers 

W(x) = lim V ^ — {- exp(f3nVx) Z(n) (8) 

n=0 

Of course, in the case Z(n) ~ exp(rt 2 ) at large n, (which is the case of REM) the above series is not that innocent. 
Here in accordance with the troubles conservation law instead of the problem of analytic continuation to non-integer 
n's we are facing formally divergent series. Nevertheless, it can be shown that such sign alternating series can be 
easily regularized, similarly to formally divergent sign alternating series 'J2'kL (—^) k a k = (l+a)( _1 * ) which at \a\ > 1 
is well defined as the analytic continuation from the region \a\ < 1. 

III. RANDOM ENERGY MODEL 
A. Definition 



The Random Energy Model is defined as a set of M = 2 states, characterized by random energies {Ei} (i — 
1,2,..., M) which are considered as independent quenched random parameters described by the Gaussian distribution 



M 



P[E 1 ,E 2 ,....,Em] = n 

1=1 

Correspondingly, the partition function of REM is 



V2^N V 2NJ 



(9) 



M 

Z = Y, ex P (-/3£0 (10) 
i=i 

which is a random quantity depending on M random parameters E\, E2, Em- The choice for the value of M — 2 N 
is motivated by the fact that the free energy of this system (as will be shown below) is extensive in lnM oc N, and 
thus it is the parameter N which plays the role of the effective "size" of the system (which is taken to infinity in the 
thermodynamic limit). The particular form 2 N (instead of say, exp(TV)) is motivated by the idea to imitate a random 
Ising system consisting of N spins (having 2^ energy states, which, of course, are not independent). 

B. "Naive" solution 

Naively, one could propose very simple derivation for the average value of the free energy of this system. Since 
the partition function, eq. (|10[) , is given by the sum of large number M of independent random terms, it could be 
approximated as follows: 



Z ~ M X exp(-^£) (11) 
Performing simple Gaussian averaging and substituting M = 2 N we get 

Z ~ 2 N exp(i7V/3 2 ) = exp(-/3iV/(/3)) (12) 

where 

fW) = -~/J-iln2 (13) 

is the free energy density of the system. Correspondingly, for the entropy we get 

S(T) - /3 2 ^/(/3) = ~\(3 2 +\n2 (14) 

Since we are dealing with the discrete system, one can immediately note that something is very wrong, as the entropy 
becomes negativ e for j3 > \/2 In 2. In fact, it turns out that there is the phase transition in the considered system 
at f3 c — \J2 In 2, such that at j3 > /3 C (in the low temperature phase) the system occupies only finite number of the 
lowest energy states. For that reason the original hypothesis of the above derivation, that partition function, eq. flTUl) . 
contains macroscopic number of random terms, turns out to be wrong at low enough temperatures. 
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C. Rigorous solution 

The result of the rigorous (non-replica) derivation of the average free energy density of REM [|| is in the following: 

-i/3-iln2, at P < /3 C = V2h[2 



fW) 



(15) 



-V21n2 , at P > 
Correspondingly, for the entropy density one gets: 



S(J3) = 



\P 2 + ln2 , at P < /3 C 







at p > p c 



(16) 



The above result for the entropy demonstrates that indeed in the low temperature phase the system effectively occupies 
only finite number of the lowest energy states. 



D. Replica approach 

In terms of the replica approach for the n-th momentum of the partition function, eq. (|10p , we get 



Z(n) = 



M 



eM~PE t ) 

n 

E 



E 



mi,...mM— 



m 1 \...m M l - 



exp 



M 



Ei-m, 



M 

^(E m *' n 



mi,...niM—0 



mi\...m,M ! 



M M 

expQiV^^m^^X: 



mi, n 



(17) 



where S(p,qj is the kronecker symbol. Reorganizing the terms, we can also get the expression which would contain 
summations over only non-zero values of ro's: 



M 



Z{n) = £ 



Ml 



E 



mi\...m k \ 



exp(±Np 2 Y, m l) «^(E m - 



(18) 



*-"k\(M-k)\ 

k=l mi,...,mfe-l a=l a — 1 

We see that at large n, 

Z{n^> 1) - exp[Cn 2 ] (19) 
which indicates that analytic continuation of the replica partition function Z(n) for non- integer n is ambiguous [5j. 



E. "RSB-magic" solution 

Heuristic RSB procedure for computing the thermodynamic limit average free energy using the above expression 
for the replica partition function, eq. (|18|) . is in the following. 

(1) All Tfi a 's are taken to be equal: mi — ... = m& = m. Correspondingly, the constrain Yla=i TO « = n turns into 
km = n which fixes k = n/m In this way the expression (|18p will contain no summations any more. Since k <C M 
we can estimate M\/(M — k)\ ~ M k — exp(/cA r ln2). Then, neglecting all pre-exponential factors (which are not 
extensive in N) the expression in eq. (fT8|) is estimated as follows: 



Z(n) - exp(-Np 2 nm + —\n2] = cxp( -PNn f(P,m) 
\2 m J \ 



where 



f(p,m) = -I^to- J-1H2 
2 pm 



(20) 



(21) 
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is a kind of the "trial free energy". The expression for f(/3,m) explicitly depends on m which is considered as an 
optimization parameter. 

(2) According to its definition the parameter m is constrained by the condition f < m < n (the fact that it must be 
multiple of n is ignored). In fact, this condition is reformulated in somewhat different way, namely: the value of the 
parameter m must be between 1 and n. Thus, after taking the limit n — > the above constrain turns into < m < 1, 
which implies that m turns into a continuous real parameter. 

(3) The parameter m is defined by the condition that the "free energy" /(/?, m), ea.(|21|). has a maximum at the 
interval < m < 1. Thus, by definition, the physical free energy of the system is 

/(/3) = max f((3,m) (22) 

0<m<l 

One can easily check that at ft < f3 c = \J2 In 2 the maximum of the function /(/3,m) in the interval < m < 1 
is achieved at m = 1, which yields f(/3) = — 1/3 — 4ln2. On the other hand, at /3 > j3 c the maximum of the 

function /(/3,m) is achieved at m — (3 c /f3 < 1, which yields /(/3) = — V21n2. Thus we see that the above "magic" 
manipulations nicely reproduce the rigorous solution, ea.([T5|). All that is quite impressive, but, of course, the above 
"derivation" can not be considered as grounded. 

F. Rigorous replica solution 

In fact rigorous replica calculations of the probability function W(x), eq.®, are rather simple. Substituting the 
expression for the replica partition function, eg. (1171) into eq.(|H]) we get 

M M 



n— mi....mM-0 ' i— 1 i— 1 

where M = 2 N and due to the presence of the kronecker symbol the summations over m^'s can be extended to infinity. 
Summing over n we can lift the constrain n = m, which provide independent summations over m^'s: 



W(af) = ^lim ( ^ V — ^ exp -iV^ 2 m? + A^zm, ) (24) 
Elementary summations of these series yields: 



i^l \mi— 



2 J ' 

VK(a;) = Jim G(N,x) (25) 



where 

r A~ / ,+0 ° . / I 



G(A, x) = \l — J di exp --AC 2 - exp(/3A( 2 ; - £)) (26) 



Taking the limit A — > ex) one easily obtain the following result (see Appendix): 



(1) In the region /3 < y/2 In 2 = /3 C 



1 , for x < -(§£+ihi2 
in,) = < (27; 
, for x > iln2 



(2) In the region (3 > V2 In 2: 



1 , for x < -V21ri2 

ir(.r) = { (28) 
, for x > -V21n2 
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According to the definition, eq.©, the above result means that in the thermodynamic limit the free energy distribution 
function of the considered model is the delta- function (which means that the system is selfaveraging) : 



In its traditional formulation the replica method procedure is in the following: first, one has to calculate the 
disorder average of the integer n-th power of the partition function, Z n — Z(n); second, one has to perform an 
analytic continuation of this function for arbitrary real or complex values of n; and third, one has to take the limit 
n — \ (to get the average free energy) or integrate over complex n (to derive the free energy distribution function). 
The third step is usually accompanied by taking the thermodynamic limit, which assumes that the system size L is 
taken to infinity. The prescription of the replica method indicates that the two limits, n — > and L — > oo, has to be 
taken simultaneously such that the product nL u (where an exponent u> defines the scaling of the free energy with the 
system size) is kept finite. 

In fact, the whole experience of the replica calculations in disordered systems shows that except for trivial cases this 
program, as it is formulated above, is never followed (for more detailed discussion of this issue see [l3j). The typical 
illustration is provided by the studies of the Sherrington-Kirkpatrick model of spin glasses [l2[ • The replica solution 
of this model which is generally believed to be correct is derived in terms of the RSB technique [l| in which all the 
above three steps, (computing Z(n), analytic continuation in n and the limits n — > and L — > oo) are performed 
simultaneously. 

In this paper an alternative replica technique has been discussed. In terms of this approach no analytic continuation 
for non-integer values of the replica parameter n is required, and instead the summation over all positive integer 
momenta of the partition function has to be performed. Earlier this method has been successfully applied for solving 
the one-dimensional directed polymer problem [t| [l(| In this paper it has been demonstrated that in terms of this 
technique the rigorous solution of the Random Energy Model takes just a few lines. Of course, real challenge would 
be to find an alternative solution for the SK model. Unfortunately, here the situation is much more complicated, as 
already at the stage of calculation of the replica partition function Z(n) (for integer n) the saddle point approximation 
is required which doesn't seem to be legitimate in terms of the present technique. In any case further systematic 
studies of the considered approach is required. 



P(f) = S(f /(/?)) 



(29) 



where the free energy f(f3) coincides with the one obtained in the rigorous (non-replica) solution, eg. (1151) . 



IV. CONCLUSIONS 
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Appendix 



Let us study the properties of the function, eq. (f2"o]) 




(A.l) 



in the limit of large N. First of all one can easily see that at x > and N ^ 1 




(A.2) 



Substituting this into eq. ([25]) . we find that 



W(x > 0) = 



(A.3) 
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At x < we can represent the function, eq. (|A.ip . as the sum of two contributions: 

G(N,x) = G 1 (N,x)+G 2 (N,x) 

where 



Gi(N,x) 



dt; ew(-±N?-ext>(-PN(\x\+tj) 



and 



G 2 (N,x) = ^—J d£ exp( --A£ 2 -exp(-/3A(|x|+0) 



Simple analysis shows that 



Gi(N, x) 



exp 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



The function G2(N,x) in the limit of large N can be estimated as follows: 



G 2 (N, x) 



+ CO 



l-exp(-(3N(\x\+0) 



rw r + °° / i \ 

1 - V 2^ 7 ^ exp(--iVe 2 - /3A£ - /3iVk|J 



-|*| 

The function <p(£) — — ^A£ 2 — /?£ has the maximum at = —j3. Thus 



(A., 



+ 00 



d£ exp(--N£ 2 - - 0N\x\ 



exp 



exp 



(-\Nx 2 ) 



for |x| < P 



(A.9) 



-/3N(\x\ - \0) , for |z| > /3 



The results, eqs. (|A.7[> . (IA.8I) and (|A.9[) . demonstrate that at N ^> 1 and negative x the function G(N,x), eq. flA.lD . 
takes the form 



G(N,x) ~ l-. 9 (A,x) 



where 



7(tf,a:) 



exp -iiVx" 4 , 



for \x\ < P 



(A.10) 



(A.11) 



exp 



-pN{\x\ - i/3) , for |s| > 



Substituting eas. (|A.10p and (jA.111) into eg. ([25]) we get 



W(a;) = lim exp —ip(N,x) 



(A.12) 



where 



^(N,x)=g(N,x)2 

Simple analysis of this expression yields: 
(a) at \x\ < P, 



N 



exp (-± Ax 2 + AH n 2) , 

p N (\x\-\p-± 



exp 



In 2 



for \x\ < P 
for |af| > P 



(A.13) 
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(b) at |x| > P, 



N 



i , for \x\ > y/2 In 2 

lim tp(N,x) = ( (A.14) 

1 -cm , for Id < V / 2hi2 



0, for \x\ > \B + ^ln2 
lim ^(i\T,a?) = { ' (A.15) 

hoo , for |ar| < \fi + ^ln2 



Substituting eqs. (|A~T4|) and (fA~T5|) into eq. (|A~T2)) we find: 

(a) at \x\ < (3, 



1 , for \x\ > V21n2 

W{x) = { (A.16) 

, for Id < V2h[2 



(b) at \x\ > (3, 



1, for \x\ > i/3+ iln2 
W(a;) = <( " (A.17) 

0, for |a;| < i/3+ iln2 

One can easily see that the above results, eq. (|A.16|l (valid for |a;| < /3) , and eq. (|A.17|l (valid for |x| > 0) are 
equivalent to eqs. (|2T|) -([2"g |) . 
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